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Abstract 

We formulate and analyze a hybrid system model that involves Volterra integral op- 
erators with multiple integrals and two types of impulsive terms. We give a constructive 
proof, via an iteration method, of existence and uniqueness of solutions. 

I. Introduction 



The class of Volterra integral equations has been traditionally used to model the 
behavior of systems with memory. The question of existence and uniqueness of solu- 
tions of Volterra integral equations is carried out by applying appropriate fixed point 
theorems; in addition to the classical works [4,8], we mention the papers [2,9, 10] that 
contain techniques related to the material of the present paper. 

In this paper, we formulate and analyze a novel class of equations that contain some of 
the characteristics of Volterra integral equations, but also fall into the general category 
of hybrid systems theory. 

The discipline of hybrid systems aims to formulate and analyze models of systems 
that combine continuous and discrete features, as well as to solve associated problems 
of control theory and game theory. In the area of ordinary differential equations, the 
standard model of hybrid systems is the class of impulsive differential equations [1] . The 
topic of Volterra integral equations with impulses has received considerable attention 
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in the recent research hterature[5, 6, 7, 11]. To the best of our knowledge, the paper[3] 
is the first to contain a general Volterra type term describing the impulses; previous 
papers contained particular types of impulses. 

In this paper, we have set out to explore the inclusion of multiple- integral Volterra 
terms and multiple types of discrete impulsive terms, with two qualitatively different 
types of impulses; the various ingredients can be combined in several different ways. It 
is expected that the present work will form the basis for further developments in the 
area of systems described by Volterra equations with multiple integral terms and also 
multiple integral-sum terms. In order to make these concepts clear at this stage, we 
shall describe the various types of integral equations that can arise as extensions and 
generalizations of the classical Volterra model. 

The standard (nonlinear) Volterra equation of the second kind for a scalar-valued 
unknown function y{t) is 

(1-1) y{t) = yo{t)+ [ f{t,s,y{s))ds 

Jo 

It is well known that general smooth phenomena with memory are not always modelled 
by single integral operators of the type /(t, s, y{s))ds, but rather they involve series 
of multiple integrals. The multiple integral series analogue of (1.1) is 

yit) = yoit) + y2-i / / fn{t,Si,S2, ,Sn, 

y(si), y(s2), • • • , y{sn))dsids2 ■•■dsn 

and a continuous solution is sought over < t < T . 

If the functions are continuous in all their arguments, satisfy bounds 

|/(t, Si,-- - ,Sn,Xi,X2,- ■ ■ ,Xn)\ < 

and Lipschitz conditions 

n 

\fnit-i Sl-i ■ ■ ■ , Sm Xi, X2, ■ ■ ■ , Xfi) fn{t-i "^1, " " " , ^n? ^1) ^2) ' ' ' ) Cn)| ^ -^'n ^ ^ \Xi 

i=l 

with the series ^n^, J2'^=o ^n+i 7^ having radius of convergence R > T, then the 

existence and uniqueness of a solution of (1.2) can be shown, in a constructive manner, 
by showing that the operator S, defined by 

°° 1 /-t rt 

o^ iSx){t) :^yo{t) + J2z;] •••/ / fn{t,Si,S2, ,Sn, 

{1.6) ri\ Jq Jq Jq 

x{si),x{s2), • • • , x{Sn))dsidS2 ■ ■ ■ dSn 
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(1.2) 



is a contraction on C(0, T; i?)(the space of continuous real- valued functions on [0, T]) in 
the norm 

(1.4) ||x|| := max e'^'^lxit)]. 

0<t<T 

If fjL is sufficiently large. The proof of the contraction property relies on the estimate 

<e"''*y'— / •••/ / \fn{t,Si,S2,- ■ ■ ,Sn,x{si),x{s2),- ■ ■ ,x{Sn)) 

^1 nl Jo Jo Jo 

-/n(i, Sl,S2,-- - ,Sn,^(si),C(s2),-- - , C{Sn))\dSidS2 ' ' ' dSn 
<x> t rt pt pt n 



n=l Jo Jo ^.^1 

oo -1 /-t pt ft ^ 

(1-5) ^'^''Y.-^ ■■■ / T.^''^ LnWx - ^w, 

n=l Jo Jo Jo ^.^1 

2^— ^nXnllx-eiU 



1 _ e-/^* ^ t 



5Z -rWilk-eil/ 



u — ' ml 

m=0 



< — j-^m+l|k-ClU 

u ml 

m=0 



and the contraction property follows from the convergence of X^^^q "^-^m+i ^ind the 
fact that limp ,00 """"^^ = 0. 

It is plain that we may assume, without loss of generality, that each function 
is symmetric with respect to permutations of the symbols (sj, Xi), i = 1, 2, ■ • • , n. i.e. 
/n(^, S(t(i), ■ • ■ ,Sa(n),a;^(i),- ■ " ,Xa{n)) = fn{t, Si, ■■■ ,Sn,xi,--- , Xn) for cvcry permuta- 
tion cr of (1,2, otherwise we could replace each by its symmetrization 

(1-6) fnit-i Si) • • • ) Sni a^i, • • • , Xfj,) •= ^ ] fn{t, 'Scr(l)) ' " " j '5CT(n)) ^cr(l)5 ' ' ' > ^(T(n))) 



where Hn permutations of (1, 2, n) , and we would have 

/ ■■■/ / fn{t,Si,S2,--- ,Sn,x{si),x{s2),--- ■,x{Sn))dsidS2---dSr 

Jo Jo Jo 



(1.7) 



f-ff 

Jo Jo Jo 



fnit, Si, S2, • • • , Sn, a;(si) , a;(s2) , • • • , a;(sn))o?sio?S2 • • • dsn- 
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Under the condition of symmetry of each (1.2) can be written as 

°° rt rSn /•S2 

yit) = yoit) + ^ / ■■■/ fnit,Si,S2r-- ,Sn, 

^—lJsn=OJsn-l=0 Jsi=0 

y{si), y{s2), ■■■ , y{Sn))dsi ■ ■ ■ dSn-ldSn- 

Another extension of the standard Volterra equation (1.1) has been introduced in [3]. 
This extension involves impulsive terms, and the impulsive Volterra equation becomes 

(1.9) y{t) = yo{t)+ I f{t,s,y{s))ds+ V g{t,Ti,y{Ti)). 

•^0 i:r,<t 

(The problems in [3] included controlled equations of the type (1.9).) 

The time-instants in (1.9) are the impulsive times. The unknown function y will 
have jump discontinuities at these instants; however, the effect of the jumps in y are 
more general than in the case of impulsive ordinary differential equations. A crucial 
feature of an impulsive Volterra equation of the type (1.9) is that a constructive proof 
of existence and uniqueness of solutions requires the use of a two-dimensional vector- 
valued metric and a concept of contractions with a 2 x 2 matrix, instead of a scalar, 
contraction "coefficient". We refer to [3] for the details. 

In this paper, we have further explored the interrelationships between continuous 
and discrete components of systems with memory. To this effect, we have considered 
Volterra equations with integral terms up to order two, discrete terms up to order two, 
and mixed discrete and continuous terms up to order two, with two kinds of impulsive 
effects, impulses of order 1, at impulsive times denoted by Tj , and impulses of order 
2, at variable times denoted by <Tj(i). The detailed formulation of this model, together 
with explanations of the terminology of higher-order impulses, is presented in the next 
section. 

The constructive proof of existence and uniqueness of solutions for our second-order 
full impulsive model requires the use of a three-dimensional vector-valued metric, and 
proof of contraction property with the role of contraction coefficient played by a matrix. 



(1.8) 



II. Properties and hypotheses 



In this section, we give some properties and hypotheses for the functions. 
Now we consider the following basic Volterra integral equation 

x{t) =xo{t) + / fi{t, x{s))ds + / / f2{t,s,si,x{s),x{si))dsids 
Jo Jo Jo 

i-1 

+ <^i(^'^i'a;(r^)) + J^G'2(t,Ti,Tj-,a;(rr),a;(rr)) 

i:Ti<t i:Ti<t j=l 
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+ J2 J] G'3(i,a,(t),T,-,a;(ai(t)-),a;(r^. ))• 

i:(Ji{t)<tj:Tj<t 

Let pij, i = 1, 2, • • • , N(j, J = 1, 2, • • • , A/o.,i be the solutions of equation 

(2.2) t = ai{t) 

Define a — {crj(-) : z = 1, 2, • • • , N^r}, and we assume that, for every z = 1, 2, • • • , N^r- 
equation (2.2) has a finite number of solutions. Let 

(2.3) p := {pij :l<j< N^,i, l<i< N^}. 
We set 

(2.4) r := {to,ti,--- ,TN^} 
and 

(2.5) I:=rUp 

The space C(0, T; R; I) is defined as the space of real- valued functions x(-) which are 
bounded and continuous on every open interval (a, /3), with endpoints a and /3 in I, and 

has limits x{a'^) := lim^ x{t), x{a~) :— lim^ >q,- x{t) at every point a e ln(0,T), 

as well as limits a;(0"'") := lim^ ,0+ x{t)i x{T~) := lim^ ,t- x{t). 

Let x{t) = ^(t), x{t~) = rji and x{ai{t)~) = /3i{t). Assume that the functions 
fii f2iGi,G2, Gs and g satisfy the following conditions: 

(HI) /i(^, s, x{s)) is continuous for < s < x e C(0, T; R; I), there exist Li > such 
that 

\fl{t,S,Xl) - fl{t,S,X2)\ < Li\xi -X2\. 

(H2) f2{t, s, si, x,y) is continuous for < si < s < t, x,y E C{0,T; R;T), there exist 
-£'21, -^22 > such that 

\f2{t,s,si,xi,yi) - f2{t,s,si,X2,y2)\ < L2i\xi -X2\ +L22\yi -y2\- 

(H3) Gi{t,Ti,x{T^)) is continuous impulse function for < Ti < t < T, x G R^, there 
exist L(3^ > such that 

\Giit,ri,r]l) - Gi{t,Ti,ril)\ < Lg^\vI -Vi\- 
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(H4) G2{t^Ti,Tj,x{T^ ),x{tj )) is continuous impulse function for < tj < Ti <t <T, 
X E R^, there exist Lg^^iLq^^ > sxidi that 

\G2{t,Ti,Tj,r]},r]j) - G2{t,Ti,Tj,T]f,r]'-) \ < Lq^iIv} - Vi\ + LG22\Vj 

(H5) G3{t, ai{t), Tj, x{ai{t)~)j x{t^)) is continuous impulse function for < Tj- < ai{t) < 
t < T, X G R"^, there exist Lg^j^, Lq^^ > such that 

(H6) g{t,s,ai{s),Tj,x{s),x{ai{s)~),x{T~)) is continuous function for < Tj < ai{s) < 
s < t < T, X E C([0, T] : r), there exist Lg^ , Lg^ , Lg^ > such that 

\g{t,s,ai{s),Tj,eit),Pl{t),r]])-\g{t,s,a,{s),r„eit),Pf{t),v])\ 

(H7) There is a positive number h such that — Tj-i > h for all z = 1,2, • • • , AT^; 

~ ^j-iit) ^ ^ ^■ll i = 2, ■ ■ ■ , A^cr, and for all t e [0, T]; whenever crj(s) < r^, 
for some j e {1, 2, • • • , N^r} and i e {1, 2, • • • , N^-}, then Tj — crj{s) > h. 

III. Existence and Uniqueness of solution of the state equation 

In this section, we will show the existence and uniqueness of solution for the nonlinear 
Volterra integral equation. 

For each collection of impulse times, we seek a solution (2.1) in the space C(0, T; R; I) 
of real valued function x{t) that are bounded and continuous on every open interval 

(a,/?), with endpoints a and /3 in I and have limits x{a'^) := lim^ x{t), x{a~) := 

lim^ x{t) at every point ct e I fl (0,r), as well as limits x{0~^) := lim^ >o+ x{t)j 

x{T-) := limt_T- x{t). 

Theorem 3.1. Suppose that hypotheses {HI) ~ {HI) are satisfied. Then for every 
r and every cr, (2.1) has a unique solution x{-) in the space C(0, T; i?; I). 

Proof. We observe that the hybrid Volterra equation (2.1) implies the following 
impulsive conditions at the times and pkf 

fc-i 

^{'^k) =^K") + Tfc, x{t^)) + G2{Tk, Tk, Tj, x{t^), x{t-)) 

i=i 

/■Tfc 

(3.1) + / ^ g{Tk,s,ai{s),Tk,x{s),x{ai{s)~),x{T^))ds 

^ i:ai{s)<Tk 

+ Y G3{Tk,ai{Tk),Tk,x{ai{Tk)~),x{T,;)) 
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(3.2) 



j-Tj<Pkl 



The interval [0,T] can be expressed as [0,T] = Uo<i<M[cti, C(i+i], for some posi- 
tive integer M, with each ai e I, and all the corresponding open intervals satisfying 
{ai, ai+i) n (ttfc, ttfc+i) = for A; 7^ /, {ai, ai+i) n I = 0. 

The solution of (2.1) can be obtained inductively as follows: equation (2.1) can be 
solved on [0, ai) as a Volterra integral equation (including double integral terms, for 
which existence and uniqueness is provided by the argument in the introduction), and 
therefore a;(aj~) can be determined. Then x{a'^) can be found by the impulsive condi- 
tions (3.1), (3.2) above. The equation (2.1) becomes a Volterra integral equation(with 
multiple integral terms, but no impulses) for the restriction of x{t) to the interval 
(q;i,q:2), and consequently x{t) can be uniquely determined over the interval [0,0:2). 
Inductively, if x{t) has been determined over the interval [0,ai), then x{a'^) can be 
determined, and x(q;^) can be found from the impulsive conditions (3.1), (3.2), so that 
(2.1) becomes a Volterra integral equation(with multiple integral terms, but no im- 
pulses) for the restriction of x{t) to the interval {ai,ai+i), which is uniquely solvable 
on {ai,ai+i), therefore x{t) will be known over the interval [0,q:;-|-i). This completes 
the induction. 

We consider the space V = C(0, T; R; I) x R^^ x C(0, T; R; I) with a vector valued 
norm defined for each {^,r)i,l3i) in V denoted by ^{t) — x{t), rji = x{Tf), (3i{t) = 
x{ai{t)-) by 

m,= sup e-^'m\ 

0<t<T 

(3.3) II^IU = ^^^e-'^^^lryil, = [?7i : < i < TV^] 

= max sup e-'^'^*W|A(t)|, P = [f3j : < j < N^] 

l<i<Aro<t<T 

We define an operator 5 on F by 

(3.4) Si^,ri,p)^ [SdU,V,m 

where : V — > C(0, T; R;I), Sa : V — > R^^ and Sm ■■ V — > C(0, T; R; I); 

= xo{t)+ fi{t,s,^{s))ds+ / f2it,s,si,C{s),^{si))dsids 
Jo Jo Jo 

i:Ti<t i:Ti<t j=l 
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(3.5) 



(3.6) 



+ / XI X] g{t,s,ai{s),Tj,^{s),Pi{s),r]j)ds 

i:ai{s)<tj:Tj<t 

+ J2 G3{t,ai{t),Tj,pi{t),r]j), 



i:ai(t)<tj:Tj<t 



= xo{ti)+ fi{Ti,s,^{s))ds+ / f2{Ti,s,si,^{s),C{si))dsids 
Jo Jo Jo 

+ ^ Gi{Ti,Ti, Vi)+YY G2{ri, Ti, rj,r]i, rij) 

i<l i<l j=l 

+ / X] Y9{ri,s,ai{s),Tj,^{s),Pi{s),rjj)ds 

i:ai{s)<Ti j=l 

l-l 

+ X YG3{Tl,a^{Tl),Tj,(3i{Tl),r]j), 



= xo{ap{t))+ fi{ap{t),s,C{s))ds 
Jo 

+ / /2(0-p(*),S,Si,C(s),C(si))rfSids 

Jo Jo 



i-1 



i:Ti<ap{t) i:Ti<ap{t) 
nCTp{t) 

+ / z2 z2 9Mt),s,ai{s),Tj,^{s),(3i{s),r]j)ds 

° i:ai{s)<ap{t) j:Tj<ap{t) 

+ X X G3{(Jp{t),ai{ap{t)),Tj,f3i{ap{t)),r]j). 

i:a-i{ap{t))<ap{t) j:Tj<ap{t) 

We shall call that Sc is continuous component of S, Sd is discrete component of 5* 
and Sm is mixed component of 5". 

Lemma 3.2. The solution of equation (2.1) is equivalent to the problem of finding 
a fixed point of the operator 5" defined in (3.5) ~ (3.7) above. 

Proof. It is clear that, if x{t) is a solution of (2.1) in the space C(0, T; R; I), then 
it is plain that the triple (^*, 77*, /?*), defined by 

eit):=xit), ry*:=x(r-), /?*(t) := x(a,(t)-) 
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is a fixed point of the operator S in the space V. 

Conversely, suppose that {$,* ,r]* , /3*) G V is a fixed point of S. Let the intervals 
[ai, ctz+i] be as in the proof of Theorem 3.1. It follows from (3.5) that ^*{t) is a solution 
of (2.1) over the time-interval [0, cti). If ai G r, it follows from (3.6) that r]l — ^*{a^), 
and then it follows from (3.5) that Cic^i) satisfies the impulsive condition (3.1) at 
t = CKi. If «! G p, it follows from (3.7) that Pi{ai) = ^*{ai{ai)~), and then it follows 
from (3.5) that ^*{<7i{ai)~^) satisfies the impulsive condition (3.2) at t = cci. Inductively, 
if ^*(t) solves (2.1) over [0, ai) and ^*{a'^) satisfies the appropriate impulsive condition, 
(3.1) or (3.2), ai t — ai, then we shall show that ^*{t) also solves (2.1) over [0,a/_|_i) 
and ^*{a'i~_^^) satisfies the appropriate impulsive condition, (3.1) or (3.2), at t = Q!/+i. 
Suppose that the points ccq, cki, • • • ,ai correspond to {tx : < A < Ao}U{cr^ : 1 < < 
IIq}. Then it follows from (3.5) that ^*{t) solves (2.1) over [0,Q;i_|_i), since the equation 
obtained from (2.1) over the interval [a;, ai-^-i) as in the proof of Theorem 3.1, contains 
no impulses in {ai, a^+i). If a^+i G r, then it follows from (3.6) that ^*(a;'^]^) = Vx^+u 
and then it follows from (3.5) that, in case a^+i < T, ^*(q;^j^) satisfies the impulsive 
condition (3.1) at t = ai+i — txq+i- If (Xi+i G p, then it follows from (3.6) that 
C*(<^/+i) ~ f^ji^ji'^i+i)~) some j > //o, and then it follows from (3.5) that, in case 
a;_i_i < T, ^*(a^t|_i) satisfies the impulsive condition (3.2) at t = pjk = q;/+i for some k. 
The inclusion is complete. 

Lemma 3.3. If > is sufficiently large, then for any (C^, ry^, /5^), (^^, r/^, /3^) of 
V, we have 

{i) \\Sc{e,v\P')-Scie,v',f3')\\^ 

< aiilie' - e'lU + «i2||^' - ^'lU + «i3||/9' - 



< aailie' - e'lU + «22||^' - r/'lU + «23||/?' - /3'lU, 

(m) \\SU^\v\P')-Smi^^v^f3^)\\^ 

< asiWe - e'lU + «32|h' - r/'lU + «33||/?' - /?'IU, 

where all constants a^- are nonnegative elements of the matrix A := {0'ij)i<i,j<3 

Proof. We will show that the operator S satisfies the contractive condition. It is 

enough to have lim^ >oo -D = 0, lim^ >oo T = and lim^ »oo S — for D — det{A), 

T = tr{A) and S = 011022 + 022033 + 033O11 - 013O31 - 023O32 - 012O21 (see Appendix 

Let L2 = max{iv2i,L22}, Lq^ = max{LG2i, Lq^J, Lg = max{Lg^, Lg^, Lg^} and 
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Lq^ = max{L(33^, Logj}. Then we have 



'y^\fi{t,s,eis))-fi{t,s,e{s))\ds 
+ f r \f2it,s,s,,eis).eisi))- f2{t,s,s,,e{s),e{si))\ds,ds 

Jq Jo 

+ \Gi{t,Ti,ril)-G,{t,n,r]^)\ 

i:Ti<t 

i-1 
i:Ti<t j=l 

+ / E E \g{t,s,ai{s),Tj,e{s),f3l{s),v]) 

° i:cri{s)<tj:rj<t 

- g{t, s, ai{s),Tj, pf{s), r]'^)\ds 

+ E E \G3{t,ai{t),Tj,Pl{t),ri])-Gs{t,ai{t),Tj,pf{t),ri^)\\. 

i:ai{t)<tj:Tj<t 

Hence, by the assumption (HI) ~ 

\\Sc{e,v\P')-Se{e,v',P')\\, 

( 1 _ p-M* + fp~i^^ 1 — p~^^^ "1 

1 _ g-M(JVt-l)ft 

e'^'* - 1 

1 _ e-AfCA^t-i)/! i_e-M(^t-i)^ 



+ <j Lg, (1 + -^^j^^ ) + LG,{Nt -l + {Nt- l)(iV, - 2) + Nt{Nt - l))e- 



f 1 _ p—fJ-Nth 1 _ -ixNth ^ 

+ r r \f2{Ti,s,s„e{s),e{si))- f2{ri,s,s„e{s),e{si))\ds,ds 

Jo Jo 

+ E l^i(^Z'^i'^i) - Gi{ri,ri,r]'^)\ 

i<l 
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i-1 

+ XlXl l^2(r/, Ti, Tj, rjl, T]]) - G2{ti, Ti, Tj,r]f, ?7|)| 

i<l j=l 

i:ai(s)<Ti 3 = 1 

- gin, s,ai{s),Tj,^'^{s), (3j{s),ri^)\ds 
l-l 



+ E El^3(T^,a,(rO,r,,/3^(rO,r^})-G3(T^,a,(rO,r,,/?2(rO,r^J)||. 

i:ai{Ti)<ri j = l 



Hence, by the assumption (HI) ~ (-^7), 

( i_e-Ma-i)'i Z(Z-l)e-^'^ l-e-MM/i 
+ \Lg^ ITTh \ + z + Lg, 



, 1 _ p-M(-'Vr(-l)fe , 1 _ p-M(A^T;-l)'i , 



\h{ap{t),s,e{s)) - hMt),s,e{s))\ds 

+ / / \f2i<7p{t),s,si,e{s),e{si)) - f2Mt),s,s,,e{s),eisi))\dsrds 

Jo Jo 

+ E |GiMt),r„r7,i)-GiK(t),r„r7f)| 

i:Ti<ap{t) 

i-1 

+ E E l<^2(o-p(t), Ti, r,-, r]l,r)]) - G2{crp{t), r^, r,-, ry,^ ryj)! 

i:Ti<<jp{t) j = l 

+ y E E I^K(*),s,ai(s),r,-,e\s),Ai(s),r/j) 

^ i:ai{s)<ap{t) j:Tj<ap{t) 

- 9i(^pit),s,ai{s),Tj,^'^(s), 01 (s) , ryl ) I rfs 
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i:ai{ap{t))<ap{t) j:Tj<ap{t) 

- G3{ap{t), ai{ap{t)),Tj,/3^{ap{t)), r]f^ 

Assume that (Ti+i(t) —ai{s) > h whenever s < (Ti{t). Then, by the assumptions (HI) ~ 
(^7), 

\\sue,v\p')-Sm{e,v^p^)L 
i^i^-^- — + L,C-f^ - — ) + L,{N^^,,,r-^ — \\\e - eh 

+ 1 _ e-M + ^c?2(^a^W + 2 1 _ e-Mfe ) 

+ LG2(^^a^(t)e + 2 ) 

+ LgCTpit)N,pit) + LG,N,^(t) ^_^-^h 1 11^ t 

f 1 _ e-M(-^V^p(t)-l)'i 1 _ g-MC^cTpCt)-!)?! ^ 



2|| 



If we take a^, 1 < i, j < 3, in the above result inequahties as same as (i), (ii) and 
{Hi) of this lemma, then (an, 021, 022, ^23) ^^315 033) — (0, 0, 0, 0, 0, 0) as fj, — > 00, and 
012,013,032 remains bounded as // — > 00 ; consequently {D,T,S) — > as ^ — > 00. 
That is, for /i sufficiently large, the operator 5" will be satisfied the contractive condition. 

We consider the following iterative scheme for the solution of (2.1): (^(0), 77(0), /3(0)) 
is an arbitrary element of V; for A; = 0, 1, 2, • • • , {^{k + l),r]{k + l),/3{k + 1)) is defined 

by 

= xo{t)+ / fi{t,s,^^k){s))ds+ / / f2{t,s,si,^^k){s),^{k){si))dsids 
Jo Jo Jo 



(3.8) 



i-l 

+ E Giit,Ti,ri(^k),i) + E E^2(*,Ti,r,-,?7(fc),i,77(fc)j) 

i:Ti<t i:Ti<tj=l 
i:ai{s)<tj:rj<t 

+ E E G3{t,a^{t),Tj,p^k),i{t)^V{k),j), 

i:ai{t)<tj:Tj<t 
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V{k+i),l 



(3.9) 



Jo Jo Jo 

i-1 

+ Gi{ri, Ti, r]^k),i) + X! X! G'2(r,, Ti, tj, ri^k),i,V{k),j) 

i<l i<l 3=1 

° i:ai{s)<Ti j = l 
l-l 

+ Yl ^<^s{ri,ai{Ti),Tj,p^k),i{^i)^V{k),j), 

i:<Ji{Ti)<Ti j=l 



P{k+l),p{^) 

= xo((7p(t))+ / /i(crp(t),s,C(/c)(s))cis 

^0 

+ / /2(o-p(i),s,si,C(fc)(s),C(fe)(si))rfsids 

Jo Jo 

i-1 

(3.10) ^ Gi{ap{t),Ti,r]^k),i)+ Y ^^'2Mt)in,rj,'n^k),i,V{k),j) 

i:Ti<crp{t) i-Ti<ap{t) 3 = 1 

rupit) 

° i:(Ti(s)<crp(t) j:Tj <(jp(t) 

+ X X <^3((Tp(t),CTi(ap(t)),rj,%),i((Tp(t)),?7(fc)j). 

i:cri(crp(t))<crp(t) j:Tj<ap{t) 

Then we have: 

Theorem 3.4. Suppose that (HI) ~ (-^7) are satisfied. Then the iterative method 
defined by (3.8) ~ (3.10) above converges to (x(-), (x(rj~), a;(r^), a;(r^), ••• ,x(r^^)), 
x{cri{-)~),x{a2{-)~), ■■■ , x(a"Ar^ (■)")), as k — > oo; the convergence of ^(fc) to x{-) is 
uniform on each closed interval [a, /3] that does not contain points of I in the interior, 
in the sense that, if wc define the restrictions ^"^^ = ^[k)i a^"^(-) = x to the interval 
[a,(3] by C^^^it) = ^^k){t) for t G [a,/3), e"4(/3) = Cik){p-), x'^^(t) = x(t) for t e [a,/3), 
a;'*'^(/3) = x{(5-), then ^J"^^ — ^ x"^ uniformly on [q;,/5]. 

Proof. By Lemma 3.3, if ^ is sufficiently large, the operator is a contraction with 
respect to the vector- valued norm || ■ ||^ on V. The contraction property with respect 
to the vector valued norm || • ||^ means that, for all x,y & V, we have \\Sx — Sy\\n < 

13 



^11 ^ ~ yWni where the 3x3 real matrix A has hm^ — >oo D — 0, hm^ — >oo ^ — 0, 

hm^ >oo S = 0. Consequently, the iterates of 5", with arbitrary initial data, converges 

to the unique fixed point of S in the topology induced on V by the vector valued norm 
II • 11^; this is a well known extension of the standard Banach fixed point theorem to 
the case of a vector valued metric, and the proof proceeds as in the standard case. 
Convergence with respect to uniform convergence on each [Tj_i,ri]. The fixed point of 
S gives the solution of (2.1) by Lemma 3.2. 

Remark. It follows from (3.5) (3.7) that ?7(fc+i),z = ^(fc+i)(rf ), (3(^k+i),i = 
i{k){'^i{t)~) for all /c > 0, so that, for A; > 1, (3.5) ~ (3.7) can also be written in the 
form 

^(k+i){t) =xo{t) + fi{t,s,^^k){s))ds+ / /2(t,s,si,C(fe)(s),C(fe)(si))rfsids 
Jo Jo Jo 

i-1 

+ Gi{t,Ti,C^k)ir~))+ X1^2(t,Ti,Tj-,^(fc)(r^),^(fc)(r-)) 

i:Ti<t i:Ti<tj=l 

+ f XI XI 9{t,s,<Ji{s),Tj,^^k){s),^{k){s),C{k){r~))ds 

i:ai{s)<tj:rj<t 

+ X X G3{t,ai{t),Tj,C(k){cri{t)~),C(k){rj~))- 

i:<7i(t)<t j:Tj<t 

Of course, if ?7(o) is chosen as ?7(o),i = ^{o){r~), P(o) is chosen as /?(o),i = ^(o)((7^(t)-), 
then (2.1) holds for aU A; = 0, 1, 2, 



Appendix I. Condition for a 3 x 3 matrix to be contractive. We shall 
derive the necessary and sufficient conditions for a 3 x 3 matrix A = [ciij]i<i,j<3 to be 
contractive, in the sense that all eigenvalues of A will have modulus strictly less than 1. 
The characteristic polynomial of A is 

(A.l) /(A) = X^-TX^ + SX-D 

where T = tr{A) = an + 022 + ^33, the trace of ^, D = det(^), the determinant of A, 
and 





an 


ai2 




an 


«13 




^22 


023 


s = 






+ 






+ 








021 


022 






^33 




^32 


033 



— <3ll(i22 + Ct22033 + CtSSO'll ~ CtlSO'Sl ~ Ct23032 — Ctl2021- 



We recall, from the theory of the Routh-Hurwitz stability criteria, that a cubic poly- 
nomial ip{u!) = Pqcj^ + piuP + P2<^ + P3, with pq > 0, will have all its roots in the open 
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left-half plane (on the plane of complex numbers) if and only if the following conditions 
are satisfied: 



(A.3) pi > 0; P3> 0; piP2 > PoPs- 

The transformation A = maps Re{uj) < onto |A| < 1. By using this transfor- 
mation into (^.1), we find that A will have all its eigenvalues in the interior of the unit 
circle on the complex plane if and only if the polynomial 

g(oj) =[1-T + S- D]oj^ + [3-T-S + 3D]oj^ 
^^'^^ +[3 + T-S- 3D]uj +[1-T + S + D] 

has all its roots in the half-plane Re{uj) < 0. After some straightforward algebra, we 
find 

(A.5) pip2 - POPS = S-4:S-8D'^ + ADT. 

Consequently, assuming pq = 1 — T-\-S — D>{), the matrix A will be contractive if 
and only if 

3 - S + 3D -T 
(A.6) l + T+5 + D>0; 

8 - 4;9 - 8Z>2 + ADT > 

Further, we note that a sufficient condition for A to be contractive is that the quantities 
S, T, D should be sufficiently small. 

In case the elements of the matrix A depend on a parameter ^ (as will be the 
case in the application to the mappings defined in Section 3 of this paper, a sufficient 
condition for A to be contractive is that (S', D, T) — > (0, 0, 0) as jx — > oo. 
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